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ABSTRACT 

New 13 CO data from the BU-FCRAO Milky Way Galactic Ring Survey (GRS) are analyzed to 
understand the shape and internal motions of molecular clouds. For a sample of more than five 
hundred molecular clouds, we find that they are preferentially elongated along the Galactic plane. On 
the other hand, their spin axes are randomly oriented. We therefore conclude that the elongation is not 
supported by internal spin but by internal velocity anisotropy. It has been known that some driving 
mechanisms are necessary to sustain the supersonic velocity dispersion within molecular clouds. The 
mechanism for generating the velocity dispersion must also account for the preferred elongation. This 
excludes some driving mechanisms, such as stellar winds and supernovae, because they do not produce 
the systemic elongation along the Galactic plane. Driving energy is more likely to come from large 
scale motions, such as the Galactic rotation. 

Subject headings: ISM: clouds — ISM: kinematics and dynamics — Galaxy: disk — Galaxy: kinematics 
and dynamics 



1. INTRODUCTION 

The nature, origin, and maintenance of the supersonic 
motions or turbulence in molecular clouds remain un- 
resolved. These motions must be 'continuously' driven, 
since otherwise they decay on a cloud crossing timescale 
~ 3 x 10 6 yr. In fact, the free-fall collapse timescale of 
dense molecular clouds and presumably the timescale for 
the decay of supersonic motions are substantially shorter 
than the likely cloud lifetimes. Such driving mechanisms 
must exert influence on the shape of molecular clouds. 

There is an argument that molecular clouds are only 
converging regions i n interstellar turbulence, and thus , 
transient structures l|LarsorJll98lt iFleck fc Clarklll98lD . 
Even in this case, the shape of structures on scales simi- 
lar to the cloud sizes will reflect the mechanisms inputing 
energy at the top of the turbulent cascade. Using new 
data from the BU-FCRAO survey, we statistically inves- 
tigate the shape of molecular clouds in the inner Galaxy. 

Possible energy sources for the interstellar turbu- 
lence ha ve recently been discussed on a theoretical ba- 
sis (s ee lElmegreen fc Scalol 120041: iMac Low &: Klessenl 
12001 . The sources may include stellar feedback (i.e. 
local events) and galactic rotation (i.e. large-scale mo- 
tions). The physical mechanism converting these ener- 
gies into turbulence has not been well understood. We 
expect that it affects the morphology of the gas struc- 
ture as well as the turbulent motions. In particular, the 
shape of molecular clouds might indicate the possible en- 
ergy source and physical mechanism. We will discuss the 
origin of the velocity dispersions as well as the shape of 
molecular clouds. 
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We use 13 CO ( J — 1 — 0) line emission data from the 
BU-FCRAO Milky Way Galactic Ring Survey (GRS; see 
Simo net al.ll200ll) . The data were downloaded from the 
GRS webpage in April 2004. They cover almost the en- 
tire region in the Galactic longitude I ~ 16° to 51°, lati- 
tude b 1.0° to +0.5°, and velocity -10 to 130 km s _1 . 

The three-dimensional (I, b, v) data cube has a fully 
sampled 22" grid in space and 0.22 kms -1 grid in ve- 
locity; the resolutions are 46" (FWHM beam size) and 
0.25 kms" 1 . 

Molecular clouds are defined as topologically closed 
surfaces with the peak main-beam brightness tempera- 
ture T m b above 4 K and the boundary temperature 2 K in 
the 13 CO data cube . The lc-noise level is AT mb = 0.48 K 
ijSimon et al J 120011. and thus, emission features above 
2 K (4er) are all significant. The peak and boundary tem- 
peratures m ay be higher than those adopted in past C O 
studies (see iScoville et al.llT98ft' iSolomon et al.lll98ft1 . 

The kinematic distances are used to estimate cloud 
sizes and hence construct samples of small and large 
clouds. Distance ambiguities are solved by utilizing the 
virial- and LTE-mass relation as a distance indicator; we 
choose the distances so that the virial- and LTE-masses 
become the most similar. The relation is confirmed for 
the clouds at terminal velocities, which do not suffer from 
the distance ambiguity. The existence of the relation at 
terminal velocities justifies the assumption that virtually 
all clouds are gravitationally bound. 

Each cloud is inspected by eye in I — b, I — v, and 
b — v maps; possibly-blended clouds (about 10%) were re- 
moved from the sample (note that this operation does not 
change the following results; and the blending affected 
very few clouds as discussed in JSJ. We also remove 
small clouds with diameter D < 5 pc in order to select 
molecular clouds rather than small clumps and to have 
clouds which were well-resolved. Past CO surveys had 
similar cut-off for cl oud identification because of their 
surve y spacing (3'; IScoville et all 1198ft ISolomon etafl 
1198ft) . The diameter is calculated as D = (AlAb) 1/2 by 
using the maximum linea r extents in the I and b direc- 
tions l)Scoville et alJll98"ft^ . Incompletely sampled clouds 
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at the edge of the cube are also removed. The final cata- 
log contains 552 clouds, most of which occupy more than 
300 pixels each in the original data cube. 

The data covers Galactic longitude I = 16 — 51°. Its 
enclosed area is about 6.0 x 10 7 pc 2 , assuming 8.5 kpc 
for the Galactocentric radius to the Sun. The average 
surface number density of our 13 CO clouds is therefore 
Nmc ~ 10 kpc -2 in the inner Galaxy (Rgc ~ 4 — 8 kpc). 
The number density is riMC ~ 100 kpc" 3 assuming 
100 pc for the thickn ess of the Galactic molecular disk 
ijSanders et alJlT984f) . These are slightly smaller than 
the values of CO molecul ar clouds (15kpc~ 2 , i.e. 1427 
clouds within 9 x 10 7 pc 2 ; iScoville et al]ll987[L 

Figure ^ shows the distribution of molecular clouds 
in v — I diagram. Four thick lines indicate the Perseus 
arm, Sagittarius ar m, Scutum arm, and 4 kpc arm (see 
iSanders et al.lll985|) . The other areas correspond to in- 
tcrarm regions. These spiral arms are identified entirely 
from the distribution of HII regions observed in radio 
recombination lines, and not from the molecular gas 
distribution. Most of massive 13 CO molecular clouds 
(> 10 5 MqJ are associated with spiral arms, but many 
less-massive clouds exist in interarm regions. This in- 
dicates either that the most massive molecular clouds 
are broken up after passing spiral arms or that they 
cool down quickly so that 13 CO emission becomes unde- 
tectable. On the contrary, the less-massive clouds survive 
across the spiral arms. 

3. COMPARISON WITH 12 CO MOLECULAR CLOUDS 

Table ^ summarizes the average properties of 13 CO 
molecular clouds. We adopted two definitions of 
characteristic size and velocity. The D and AV 
are th e diameter and full v elocity width of the 
clouds l|Scoville fc Sanders! Il987|) . The 5" and ay are 
the intensity-weighte d dispersions of size and velocity 
(Solo mon et al.l 11987^ The virial mass w as calculated 



with the equation in lSolomon et alJ l|1987j) . 

The diameter D is consis tent with the CO val ues both 
on the average (18.1 pc in IScoville et al.lll987fl and on 
the m ass-weighted average C40 pc in IScoville fc Sanders! 
1987). Therefore, the extensions of 13 CO molecular 
clouds are almost the same as those of CO clouds. The 
velocity dispersion cry, however, is 1.4 time s smaller than 
the o ne in past CO studies (1.4 km s" 1 in Scovill eet al.l 
1987). The a v -D relation is a v = 0.23Z? 054 , whose co- 
efficient is ~ 1.4 times sm aller than the one i n the CO 
studies (cr v = 0.31D 55 in IScoville etld1ll987Th A natu- 
ral interpretation may be that the central dense regions 
of molecular clouds were saturated in the past CO stud- 
ies. The S, i.e. dispersion, is also s maller by a factor of 
several than the one in CO studies l|Solomon et alJll98^. 
; see their Figure 1). Further discussion of the size and 
mass distributions will be presented in Sawada et al. (in 
preparation) . 

The averaged masses are different by an order of mag- 
nitude between the mass-weighted average and normal 
average (Table . This indicates that most molecular 
gas is in giant molecular clouds (> 10 5 M Q ), however, 
the majority of molecular clouds has only ~ 10 4 M Q . 
These averaged masses are consis tent with the ones 
derived in the past CO studies ijScoville et alJ Il987t 
IScoville fc Sandersll987(K if the 1.4 times difference of a v 



is taken into account. The sound speed corresponding to 
the molecular gas temperature of 10K is ~ 0.2kms _1 . 
The much larger velocity widths (AV and cry) indicate 
that the internal gas motions are supersonic. 

4. ELONGATION 

The position angles and axis ratios of molecular clouds 
are calculated by diagonalizing the moment of inertia 
matrix as, 
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where TJy is the brightness temperature at a pixel 
and otij and are the distances from the emission cen- 
troid to the pixel in the I- and ^-directions. Re is a rota- 
tion matrix with the rotation angle 9. The axis ratio and 
position angle of clouds correspond to (I xx / 'lyy) 1 ^ 2 and 9, 
respectively. The position angle is defined from the pos- 
itive 6-direction going counterclockwise to the positive I. 
The direction along the Galactic plane is 9 = 90°. 

Figure [3 shows the histograms of axis ratio and posi- 
tion angle. The axis ratio has the peak at about 1.8, in- 
dicating that the molecular clouds are significantly elon- 
gated. The population of round clouds, having the axis 
ratio of 1, is small. The position angle peaks at 9 ~ 90° 
(i.e. along the Galactic plane). Clouds with 9 ~ 90° is 
about two times more populated than those with 9^0° 
or 180°. Therefore, molecular clouds are elongated pre- 
dominantly in the direction of the Galactic plane. We 
separate small clouds (D < 15 pc) from large clouds 
(D > 15 pc) in Figure As well as the large clouds, 
the small clouds are predominantly elongated along the 
Galactic plane. These clouds are much smaller than the 
width of the gas disk (about 100 pc) and stellar disk 
(about 600 pc). 

5. SPIN 

It is possible that the elongation is supported by the 
spin of molecular clouds. In order to evaluate this, we 
measured the velocity gradient in molecular clouds (pro- 
jection on the sky) and the position angle perpendicular 
to the maximum velocity gradient, i. e. parallel to the 
axis of possible spin (Figure |2| • The velocity gradient 
is estimated from the observed velocity field by fitting a 
plane. The position angle ip is oriented with 0° at +b, 
increasing toward +1. Typical velocity gradients dv/dr 
projected on the sky are about 0.15 kms -1 pc -1 . We 
did not correct for the apparent velocity gradient due 
to the projected LSR velocity over the diameter of the 
clouds, since this gradient is at most ~ 0.04 kms" 1 pc" 1 , 
assuming a rotation velocity of 220 km s _1 . Prograde and 
retrograde spins with respect to the Galactic rotation are 
indicated with different patterns, and are equally popu- 
lated. 

Figure |3| right can be directly compared to Figure |21 
right. If spin supports the elongation, the ip should be 
distributed similarly to the 9; their peaks should be re- 
lated as V'poak ~ # P cak — 90°. The figures, however, show 
that the spin is randomly oriented and has no dominant 
peak. Figure0]shows the difference of the position angles 
of the elongation axis and spin axis. This figure shows 
no significant peak at 90°, compared with Figure right 
and considering the statistical error ~ TV 1 / 2 . Therefore, 
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the cloud elongations are not supported by internal spins, 
but should be related to the internal velocity dispersions 
(supersonic motions). 

It is noteworthy that confusion (blending) in cloud 
identification is likely to align the apparent spin axis per- 
pendicular to an elongation - it is most likely that two 
clouds are slightly offset in space and velocity. The op- 
posing result strongly indicates that the blending is not 
important in our sample. 

6. DISCUSSION 

The supersonic motions within molecular clouds must 
be regenerated continuously; otherwise they decay 
rapidly within a cloud crossing time. Mechanisms for 
driving the supersonic motion must also account for the 
cloud elongation along the Galactic plane. In light of 
the new results, we discuss a possible energy source 
and mechanism for driving supersonic motions and cloud 
elongations. 

6.1. Energy Source 

Among several sugg ested energy sources (see 
Elm egreen fc Scald H)04), two sources, i.e. stellar 
feedback and galactic rotation, are most often consid- 
ered. Most other sources extract their intrinsic energies 
from these two. 

Stellar feedback includes protostellar outflows, stel- 
lar winds, ionizing radiation, and supernova explosions. 
These feedback mechanisms, however, release energy 
isotropically or in random directions. Such processes are 
unlikely to cause the preferential orientation of molecular 
clouds. They would cause a random orientation. 

Supernovae are the most e nergetic feedback mechanism 
l|Norman fc Ferraral 1199(1 iMac Low & Klessenl 120041) . 
There are some arguments for supernovae-induced dense 
gas formation and supersonic turbulence. If supernova 
blastwaves, however, compress surrounding material and 
induce molecular cloud formation, their round expansion 
will produce no preferential orientation. The blastwaves 
would escape from the Galactic disk predominantly in 
the direction normal to the Galactic plane. Then, most 
molecular clouds would be elongated perpendicular to 
the plane. The stellar feedback is unlikely to be the dom- 
inant cause of the elongation of molecular clouds. 

Evidence has been accumulated for the presence 
of molecular clouds w i thout star formation activ- 
ity lIMoonev & Solomon! Il98^: IScoville fc Goodl ll989T: 
iWilliams fc Blitd 11998ft . Stellar feedback cannot sup- 
ply energy to such clouds. High mass star-forming re- 
gions are localized around spiral arms in the Galaxy 
l)Downes e1^_alj[19_gp^J^iK iers et al.ll985t) and in external 
galaxies (|Scoviile et al]l200ll) . It is also unlikely that the 
energy supply from stellar feedback suffices to maintain 
supersonic motions in interarm molecular clouds. 

Since the cloud elongation has the preferential direc- 
tion, the energy source is likely to have a similar preferred 
direction. Galactic rotation has an rotation axis, and is a 
likely energy source. Galactic rotation provides an enor- 
mous reservoir of energy; it can maintain the elongations 
and supersonic motions over the Galactic age. 

Even if molecular clouds are only transient structures 
in large-scale interstellar turbulence, the anisotropic 
shapes should indicate the origin of the turbulence. It is 
still likely that the Galactic rotation is the main source 



of the turbulant energy. 

6.2. Driving Mechanism 

How might the rotation energy be transfered to 
random motion in molecular clouds? Energy trans- 
fer from rotation to random motions generally occurs 
in a differentially rotating disk, if mass elements or- 
biting at slightly differe nt radii pull on each other 
l|Sellwood fc Balbusl IT991 . iSellwood fc Balbusl l|1999fi 
pointed out that many types of viscous stress, including 
Reynolds, Maxwell, and gravitational stresses, can cause 
the energy trans fer from rotation to ran dom motion. 

For example, Uog fc Ostrikeil l)1988j) found with an 
analytic calculation that cloud-cloud gravitational in- 
teractions (gravitational stress) can transfer the or- 
bital energy of clouds into random motions between 
clouds [but they did not disc uss internal motions] . 
iWada. Meurer fc Normanl (|2002f ) used hydrodynamical 
simulations and showed that the rotational energy could 
cascade into smaller scales (down to their resolution 
limit, i.e. the size of giant molecular clouds), owing to 
local tidal force among clouds, filaments, and voids. Per- 
haps, this cascade proceeds down to the internal motions 
of the clouds. 

If tidal gravitational torques exerted by passing clouds 
is the cause of the velocity dispersions, the preferred 
elongation of molecular clouds would be produced natu- 
rally. The mean separation of molecular clouds is about 

1 /3 

V n MC ~ 200 pc. This is larger than the thickness of the 
Galactic molecular disk 100 pc. Tidal radius of molecular 
clouds should be at least a few times larger than their 
diameter, and thus is not negligible compared with the 
disk thickness. Tidal interaction should occur predomi- 
nantly in the Galactic plane; it may generate asymmetric 
velocity dispersions which cause the elongation along the 
Galactic disk. 

A slight exchange of angular momentum may occur 
during tidal interaction and spin up clouds. The spin axis 
might become perpendicular to the elongation, however, 
this effect should be little. In case that two clouds pass 
by at the distance d = 100 pc and velocity v = 6kms _1 
(typical cloud-cloud velocity dispersion), the tidal inter- 
action would increase the velocity gradient within a cloud 
by at most v/2ird ~ 0.01 kms -1 pc -1 if the two clouds 
are strongly coupled. This is negligible compared with 
the observed velocity gradient. In addition, random en- 
counters would suppress the increase of angular momen- 
tum. The spin-up is not significant even if the tidal in- 
teraction is the cause of supersonic motions. 

We discussed only gravitational interactions between 
molecular clouds. Some other density fluctuations, such 
as filamentary structure in gas, would also strengthen 
the tidal field around clouds. There are, of course, some 
other possible mechanisms introduced by magnetic field, 
surrounding diffuse gas (HI and H2), and gradients in 
stellar potential. Further discussions would require the- 
oretical researches. 

7. SUMMARY 

We showed that 13 CO molecular clouds are elongated 
predominantly along the Galactic plane. The elongations 
are not supported by the spins of molecular clouds, but 
should be supported by supersonic velocity dispersions in 
molecular clouds. In order to avoid a rapid decay of the 
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supersonic motions, there must be some mechanisms to 
drive them continuously. The driving mechanisms must 
also account for the elongations of molecular clouds. The 
driving energy is likely to come from large-scale motions, 
such as the Galactic rotation, and not from stellar feed- 
back. 
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TABLE 1 

Properties of 13 CO Molecular Cloud 



Parameter 


mass- weighted mean 


mean 


D (pc) a 


38 


15 


S (pc) b 


8 


3 


AV (kms- 1 ) 1 


11.5 


5.1 


ov (kms- 1 ) d 


2.2 


1.0 


Tpeak (K) e 


11 


8.9 


M VT (Mq)/ 


1.2 x 10 5 


1.3 x 10 4 



Note. — Mean parameters are calculated among molecular clouds with D > 5 pc 

a Diameter of molecular cloud 

b Intensity- weighted size 

c Velocity width at the 2 K level 

d Intensity- weighted velocity dispersion 

°Excitation temperature (calculated from main beam temperature and cosmic microwave background temperature (2.7 K)) 
f Virial mass 




50 45 40 35 30 25 20 15 
Galactic Longitude [deg] 

Fig. 1. — v — I diagram of the distribution of all identified molecular clouds. Thick grey lines indicate four major spiral arms, i.e. Perseus 
arm, Sagittarius arm, Scutum arms, and 4kpc arm (from left to right). These spiral arms arc identified entirely from radio recombination 
observations, and not from CO observations (Downes et al. 1980; see also Sanders et al. 1985). 
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Fig. 2. — Distribution of axis ratios (left) and position angles (right) of molecular clouds. The values are calculated using the moments 
of inertia. Position angle 9 is defined from the +6 direction going to positive I. Hence 9 = 90° is the direction along the Galactic plane. 
Large (D > 15 pc) and small (D < 15 pc) clouds are indicated with different patterns and illustrated accumulatively. 




Fig. 3. — Distribution of velocity gradients (left) and position angles of spin (right) of molecular clouds. Position angle ■)/> is defined as 
in Figure Hl(ngto), but ranges between —90° and 90°. Hence, ij> = 0° is the direction perpendicular to the Galactic plane. The right figure 
can be compared directly with Figure |5] (right). Prograde and retrograde spins are indicated with different patterns. 



